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ABSTRACT 


In  this  thesis  we  present  two  simple  models  of 
nuclear  rotations. 

The  first  model  consists  of  a  mass-point  rotating 
in  a  stationary  potential.  This  potential  was  deducted 
from  the  energy  relation  of  the  Harris  model.  The  results 
of  the  quantized  version  of  this  model  are  compared  to 
experiments . 

The  second  model  is  the  field  theoretical  equiva¬ 
lence  of  particles  confined  to  a  circle  of  a  fixed  radius. 
The  particles  interact  through  an  interaction  dependent  on 
their  relative  angle  and  their  angular  momentum. 

Approximations  of  the  energy  eigenvalues  of  the 
first  model  do  fit  quite  well  to  the  experimental  results 
for  the  ground  state  rotational  bands  of  typical  rota¬ 
tional  nuclei.  Due  to  the  approximation  method  the  fit 
is  not  so  good  for  the  lowest  excited  states  of  each 
nuclei.  The  predicted  values  of  the  moment  of  inertia 
and  the  quadrupole  moment  are  a  bit  lower  than  the  experi¬ 
mental  ones. 

The  second  model  gives  a  very  good  fit  to  experi¬ 
mental  results  both  for  high  and  low  states  of  angular 
momentum. 
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CHAPTER  1 


INTRODUCTION 

The  study  of  rotational  spectra  of  nuclei  has  given 
a  valuable  insight  into  collective  and  single  particle 
motion,  and  their  coupling  in  the  nucleus. 

Attempts  were  made  as  early  as  1938  by  Teller  and 
Wheeler  to  detect  rotational  energy  spectra  of  few  nuclei, 
but  as  they  expected  the  rotation  to  be  a  property  of  the 
whole  nucleus,  thought  of  as  a  rigid  rotor  with  the  corres 
ponding  classical  moment  of  inertia,  and  due  to  inadequate 
equipment  this  search  was  without  any  success. 

In  1951  Bohr  put  forward  the  theory  that  rota¬ 
tional  excitation  modes  of  the  nucleus  were  the  necessary 
consequence  of  a  strongly  deformed  equilibrium  state 
(spontaneous  rotational  symmetry  breaking)  (Bohr,  51) . 

The  same  nuclei  should  therefore  show  large  quadrupole 
moments . 

Shortly  after  rotational  bands  were  identified,  in 
observation  of  level  sequences  proportional  to  £(£+1) 

(£:  angular  quantum  number)  (Bohr  and  Mottelson,  75  p.  2). 

Most  of  the  early  observed  rotational  bands  were 
in  agreement  with  the  energy  expression: 

2 

E(£)  =  ft  (y  1)  (£  =  2,4,6...  even  parity) 

(I:  moment  of  inertia) 

due  to  the  fact  that  highly  excited  states  were  not 
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available  at  this  time. 

In  nearly  all  cases  the  moment  of  inertia  was 
smaller  than  predicted  by  assuming  a  rigid  nuclear 
structure  (Bohr  and  Mottelson,  75  p.  74). 

Bohr  explained  this  deviation  from  the  rigid  moment 
of  inertia  by  assuming  that  only  a  fraction  of  the  nuclear 
matter  was  participating  in  the  rotational  motion  around 
a  stationary  nuclear  core. 

Soon  it  was  discovered  that  all  nuclei  with 
150  <  A  <  186 ,  A  >  224  and  some  few  others,  had  rotational 
spectra  (Mariscott  and  Scharf f-Goldhaber,  69). 

Many  nuclei  were  found  to  violate  the  i(t+l)  depen¬ 
dence.  With  improving  methods  to  excite  high  angular 
momentum  states  in  nuclei  it  was  found  that  in  nearly  all 
cases  the  level  spacing  was  denser  than  predicted  by  the 
i(t+l)  rule. 

As  a  consequence  of  this  it  was  suggested  that  a 
term  of  the  form: 

-  B£2  (JL+1)2 

should  be  added  to  the  energy  expression  to  account  for 
rotation  vibrational  interaction  (Scharf f-Goldhaber  et.al., 
76  p.  245) ;  however  with  this  additional  term  the  energy 
expression  could  not  be  used  to  fit  experimental  results 
for  neutron  rich  nuclei,  in  which  the  rotational  and 
vibrational  level  spacing  is  approaching  the  same  order 
of  magnitude  with  increasing  N,  neutron  number  (Mariscott 
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and  Scharf f-Goldhaber ,  69  p.  1865). 

In  spite  of  nuclei  with  widely  differing  N,  z, 

* 

E(£=z)  Mailman  was  able  to  show  that  if  the  energy  ratios 
E^/E^r  Eg/E^  ^or  even-even  nuclei  were  plotted  against 
E^/E^  he  would  obtain  two  universal  curves  (Scharff- 
Goldhaber  et.al. ,  76  p.  269).  This  seems  to  indicate 
that  the  same  dynamical  features  govern  nuclei  of  all 
magnitudes  of  deformation. 

One  of  the  first  models  of  nuclear  rotation  was 
the  hydrodynamical  one  by  Bohr  and  Mottelson  ( Scharf f- 
Goldhaber  et.al. ,  76  p.  246).  It  assumed  f(£+l)  energy 
dependence  and  was  therefore  inadequate  to  describe  the 
energy  levels  for  high  Z  for  most  nuclei.  Moreover  in 
this  model  one  has  to  assume  too  high  values  for  the 
deformation  parameter  3  in  order  to  predict  values  for 
the  moment  of  inertia  divided  by  the  quadrupole  moment 
comparable  to  experiment  (Prakash  et.al.,  73  p.  2475). 

If  rotational  bands  are  assumed  to  be  described  by 

„  _  r?Ui+D 

E  _  21 - 

then  the  decrease  in  level  spacing  with  increasing  Z,  must 
be  due  to  increasing  moment  of  inertia  I.  It  has  been 
observed  that  for  most  nuclei  I  is  approaching  the  value 
for  a  rigid  nucleus  with  increasing  Z  (Mariscott  et.al., 

69  p.  1866) . 

Many  models  have  been  constructed  to  account  for 


this  increase  in  I. 


An  early  model  (3  stretching  model)  tried  has  the 


energy  expression: 


VB>  =  I  C(B£ 


£(£+1) 

2I(3£) 


I  is  here  in  units  of  'ft  ,  3  is  the  deformation  parameter. 
To  derive  the  value  of  3^ r  the  equilibrium  condition 
8E^/93^  =  0  is  used  (Mariscott  et.al.,  69  p.  1866). 

This  model  gave  a  good  fit  to  rotational  bands  of 

2 

strongly  deformed  nuclei,  if  one  assumed  1^3  (us  in  the 
hydrodynamical  model)  but  was  inaccurate  for  others 
(Mariscott  et.al.,  69  p.  1866).  One  disadvantage  of  this 
model  is  that  later  studies  of  the  nucleus  (muonic  x-rays 
and  isomer  shift  in  Mossbauer  effect)  have  shown  that  3 
does  not  increase  as  much  with  l  as  this  model  would 
indicate  (Mariscott  et.al.,  69  p.  1866). 

A  very  successful  one  has  been  the  Variable  Moment 
of  Inertia  model  (V.M.I.  model).  The  energy  according 
to  this  model  is  given  by: 


VJ) 


=  lc(I-Io)2 


1  j(ji+i) 

2  I 


with  the  equilibrium  condition: 

9E  (I )  =  n 

91  u 

Here  the  moment  of  inertia  I  is  treated  as  a  general  vari¬ 
able,  instead  of  the  deformation  parameter  (Mariscott  et. 
al.,  69  p.  1867).  This  model  has  justified  the  Mailman 
curves,  and  given  a  good  energy  fit  closer  to  the 
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spherical  region  than  the  3  stretching  model.  The  success 
of  the  V.M.I.  model  may  be  related  to  the  fact  that  noth¬ 
ing  is  stated  about  how  I  depends  on,  for  example,  the  3 
stretching  of  the  nucleus  nor  about  the  decrease  in  pair¬ 
ing  energy  with  increasing  £. 

Another  model,  the  cranking  model,  has  been  shown 
to  give  the  same  results  as  the  V.M.I.  model,  when 
expanded  to  the  fourth  power  of  the  angular  velocity  to 
by  perturbation  theory  (Mariscott  et.al.,  69  p.  1870). 

The  cranking  model  considers  a  rotating  deformed 
self-consistent  potential  well.  The  hamiltonian  in  the 
intrinsic  nuclear  frame  is 


H  “  toL 
o  z 

H  is  the  static  hamiltonian  of  the  nucleus, 
o 

4 

The  energy  expression  up  to  to  in  the  lab.  frame  is 

E  ,  =  -  to2  (I  +  3Cco2) 
rot  2  o 


(M&+1)  )**  =  co(I  +  2Co2) 

(Harris,  65).  Expression  for  E  t  independent  of  to  can  be 
obtained  from  these  two  expressions.  It  will  have  the 
same  form  as  the  energy  formula  of  the  V.M.I.  model. 

Harris  has  shown  that  the  energy  of  a  rotational 
band  can  be  approximated  by  a  self-consistency  approach 
to  any  power  of  to 

E  =i  to2  (I  +  3Cto2  +  5Dto4  +  7Fto6  +  ...) 
rot  2  o 
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U  U+l)  )  ^  =  03  (I  +  2C0)2  +  3DO)4  +  4F(06  +  ...) 

(Harris,  65  p.  B511)  .  For  the  even-even  rare  earth  nuclei 
Harris  was  able  to  fit  the  energy  levels  within  0.76%  rms 
with  just  a  two  parameter  fit,  while  the  experimental  error 
is  0.3%. 

Gupta  has  considered  an  interesting  model,  in  which 
an  axially  symmetric  prolate  deformed  nucleus  is  rotating 
around  an  axis  perpendicular  to  its  symmetry  axis.  In 
this  model  the  centrifugal  stretching  is  counteracted  by 
a  harmonic  restoring  force  (Gupta,  69). 

This  model  gives  satisfactory  fit  to  experimental 
values  of  the  rotational  energy  of  even-even  nuclei,  though 
not  as  good  as  the  Harris  model,  but  instead  it  gives  a 
suggestion  for  the  internal  mechanism  in  a  nucleus  at  the 
classical  level. 

By  increasing  the  strength  of  the  restoring  force, 
Gupta  is  able  to  reproduce  the  characteristics  of  nuclear 
vibration  of  two  different  modes  (ellipsoidal  and  axial 
symmetric) . 

This  model  was  later  renewed  by  adding  to  it  a 
rotationally  invariant  core  to  better  account  for  the 
internal  structure  of  the  nucleus  and  give  the  possibility 
that  the  moment  of  inertia  could  be  expanded  in  a  series 
in  terms  of  m  (Trainor  and  Gupta,  71)  basically  different 
from  the  expression  for  I  in  the  hydrodynamical  model. 


This  model,  named  the  governor  model,  gives  satis¬ 
factory  values  of  the  moment  of  inertia  and  the  quadru- 
pole  moment  in  contrast  to  the  hydrodynamical  model 
(Prakash  et.al.,  73). 

The  governor  model  has  been  quantized  by  Prakash 
et.al.  resulting  in  a  closed  form  expression  for  the 
energy  (Prakash  et.al.,  73)  that  gives  a  good  fit  to 
experimental  values  of  the  energy  as  before. 

It  has  been  observed  when  the  moment  of  inertia 
of  several  nuclei  is  plotted  versus  the  angular 
velocity,  that  the  resulting  smooth  curves  "bend  back" 
for  some  high  £  values,  the  graph  comes  triple  valued. 
Neither  the  V.M.I.  model  nor  the  cranking  model  according 
to  perturbation  theory  can  explain  this  phenomenon. 

To  account  for  the  back  bending  and  some  other 
"irregularities"  in  the  moment  of  inertia  a  model  that 
deals  with  the  interactions  between  the  nucleons  directly 
has  to  be  constructed  (for  back  bending  see,  for  example, 
Scharf f-Goldhaber  et.al.,  76  p.  282). 

One  model  that  deals  with  those  difficulties  is 
constructed  from  the  second  quantized  version  of  the  crank 
ing  model  hamiltonian 


H  =  H 


o 


oL 


Ho^ 


l 

ij 


t  1  v 

T  .  .  C  .  C  .  +  x  / 

13  1  3  A 


i  jk£ 


t  t 

V.  -C.'c'.C,  c0 

I3k£  1  3  k  £ 


Then  a  canonical  transformation  is  applied  to  write  H  in 

o 

terms  of  weakly  interacting  quasiparticles. 


The  inertial  forces  introduced  by  the  external 

cranking  potential ,  - wL  , become  very  strong  for  states  of 

z 

high  angular  momentum  so  perturbation  methods  are  not 
valid . 

This  model  instead  deals  with  them  exactly,  so 
shape  transitions  can  happen  at  high  %  resulting  in  back 
bending.  Moreover  this  model  is  not  limited  to  axial 
deformations  and  is  self-consistent,  as  it  relates  the 
energy  spectrum  directly  to  the  nucleon-nucleon  inter¬ 
actions  (Scharf f-Goldhaber  et.al.,  76  p.  301-312).  In 
1976  only  very  few  calculations  with  this  model  including 
realistic  forces  had  been  done,  but  the  preliminary 
results  were  promising  (Scharf f-Goldhaber  et.al.,  76, 
p.  314)  . 

Several  other  models  and  methods  have  been  tried  to 
describe  and  explain  nuclear  rotation,  but  these  examples 
served  to  show  the  mainstreams. 

In  the  first  part  of  this  thesis  we  want  to  take  a 
new  approach  to  the  problem  of  nuclear  rotation. 

It  is  known  that  the  energy  expression  of  the  Harri 

model : 

E  =  (I  +  3 Coo2 ) 

2  o 

describes  the  rotational  energy  of  even-even  nuclei  with 
good  accuracy,  as  mentioned  before. 

Is  it  possible  to  find  or  construct  a  one  particle 
potential  classically  and  quantum  mechanically  that  will 
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give  the  same  expression  for  E?  If  so,  we  can  think  of 
this  potential  as  the  average  potential  that  the  valence 
nucleons  are  subjected  to,  as  all  experiments  seem  to 
indicate  that  just  a  part  of  the  nuclear  material  is 
participating  in  the  rotation. 

This  is  done  in  the  first  part  (Chapters  2  and  3) , 
and  expectation  values  of  the  nuclear  radius  R,  moment  of 
inertia  I  and  the  quadrupole  moment  Q  are  compared  to 
experimental  values  and  other  models. 

Thieberger  (Thieberger,  70)  has  introduced  a  class¬ 
ical  model,  which  has  the  same  energy  expression  as  the 
Harris  model: 

Two  water  reservoirs  are  connected  to  each  other 
through  a  narrow  pipe.  The  smaller  one  is  rotating  around 
the  heavier  stationary  one  in  a  uniform  gravitational 
field. 

A  model  with  similar  structure  has  been  considered 
by  Bassichis  and  Foldy  (Bassichis  and  Foldy,  64).  Parti¬ 
cles  are  moving  on  a  circle,  of  a  constant  radius,  inter¬ 
acting  through  a  two  body  force  symmetric  in  the  angular 
variable  0^-0^.  If  in  the  second  quantized  version  we 
just  allow  for  the  two  lowest  energy  levels,  then  we  have 
the  lowest  level  of  zero  angular  momentum  and  a  degenerate 
level  of  opposite  and  equal  angular  momentum. 

In  the  second  part  (Chapter  4)  of  this  thesis  we 
search  for  an  interaction,  with  the  Thieberger  and  Harris 
models  in  mind,  that  would  when  introduced  to  the  second 


quantized  model  give  a  similar  energy  relation  as  the 
Harris  model. 

The  second  part  is  concluded  by  inspection  of 
several  symmetries  conserved  by  this  model. 


CHAPTER  2 


SEARCH  FOR  A  CLASSICAL  POTENTIAL  TO  EXPLAIN 
NUCLEAR  ROTATION  ENERGIES 

In  this  chapter  we  will  deduce  a  potential  that  can 
be  used  in  a  simple  classical  model  of  nuclear  rotations, 
from  expressions  that  have  been  used  successfully  to  des¬ 
cribe  rotational  energy  levels  in  many  even-even  nuclei. 

In  some  models  the  variable  moment  of  inertia  has 
been  accounted  for  by  assuming  that  the  centrifugal  force 
is  counteracted  by  a  harmonic  restoring  force,  as  for 
example  in  the  governor  model  by  Trainor  and  Gupta  (Trainor 
and  Gupta,  71) .  But  to  our  best  knowledge  nobody  has  pro¬ 
ceeded  the  same  way  we  will. 

It  has  been  shown  by  Harris  (Harris,  65)  and  others 
(Mariscott  et.al.,  69)  that  the  following  formulas  for  the 
energy  E  and  the  moment  of  inertia  I  : 

E  =  i  o2(I  +  3Co)2 )  +  E  (2.1) 

2  o  O 

I  =  (I  +  2Cco2)  (2.2) 

o 

I  :  moment  of  inertia  of  the  ground  state 
a)  :  angular  velocity 
C  :  higher  order  inertial  parameter 

do  give  an  adequate  description  of  the  rotational  energy  of 
a  wide  range  of  nuclei.  By  combining  the  two  equations 
(2.1)  and  (2.2)  we  get: 
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12 


12  4 

E  =  (I  +  CO)  )  +  E 

2  o 


(2.3) 


Now  we  want  to  find  a  model  that  has  the  following 


equation  of  motion: 


mr  -  mrt{>2  +  ~  V(r)  =  0 

o  r 


(2.4) 


V(r):  central  potential 
(J)  :  angular  variable 

which  when  in  equilibrium  (that  is  when  r,r  =  0)  will  lead 
to  the  energy  expression  (2.3). 

In  equilibrium  the  equation  of  motion  (2.4)  will 
have  the  form: 


I2  3V 

mr*  “  sr 


(2.5) 


•  2 


which  when  solved  for  <£  gives 


=  -1-  A.  v 

Y  mr  3r 


(2.6) 


The  expression  (2.2)  for  the  moment  of  inertia  will 
give  when  divided  by  the  mass  m  : 


2  ,  2  ,  2C  r 2 . 

r  =  (r  +  —  i  ) 
o  m  ^ 


(2.7) 


.  2 

If  we  use  expression  C2.6)  for  <j>  in  (.2.7)  we  get  the 
following  differential  equation  for  the  potential 


eq 

and  therefore: 
3 


r2  +  2C  ~ 
o  m2r  3r 


V  (r) 


r=r  (i) 


3  r 


V  (r) 


eq 


2 

m  r  ,  2 

-2c  (r 


2. 

ro> 


(2.8) 


r=r  ( (j> ) 


- 
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This  equation  has  the  exact  solution: 

2 

v(r)  =  ic  (r4  ~  2ror2)  ;  (2.9) 

The  system  that  has  (2.4)  as  an  equation  of  motion  has 
the  following  energy  expression  when  in  equilibrium: 

E  =  [|  mr24.2  +  V (r )  ]  .  (2.10) 

eq  2  r=r(i) 

If  we  use  equation  (2.9)  for  V (r)  in  this  expression  for 
E  we  get: 

2 

2*2  ,  m  /  4  _  2  2.  _  * 

Eeq  h  mr  ^  +  8C  r  2ror  ^  *  (2.11) 

Now  from  the  equilibrium  condition  (2.6)  we  have  the 

2  •  2 

following  relationship  between  r  and  4>  : 


•2  1  3  T7,  .  m  ,  2  2. 

*  =  E?  37  v(r)  “  2C  (r  -  ro> 


(2.12) 


therefore : 


2  2C  -2  2 

r  =  —  d>  +  r 
m  Y  o 


(2.13) 


2  •  2 

when  this  expression  of  r  in  terms  of  <J>  is  used  in  (2.11) 


we  get: 


eq 


1  2-23  *4 

r  mr  Cd) 

2  o  2 


1  4  1  .23*4 

5“rr*=yI<J>  +fC(f>  +E^  (2.14) 

8C  o  2  o  2  o 


which  has  exactly  the  form  we  wanted  (2.1) . 

So  the  conclusion  is  that,  the  potential  that 


classically  does  give  the  wanted  energy  dependence  (2.1) 

2 


:2  . 

on  (J)  is: 


tt  /  \  m  ,4  _  2  2 . 

V(r)  =  8C  (r  "2ror  > 


(2.9) 


CHAPTER  3 


QUANTIZED  MODEL  FOR  NUCLEAR  ROTATIONS 

In  this  chapter  we  will  quantize  a  simple  model 
involving  the  potential  (2.9)  which  was  shown,  in  the  last 
chapter,  to  be  able  to  give  the  successful  classical 
energy  expression  for  nuclear  rotation  (2.1). 

Then  the  energy  eigenvalues  of  the  quantized  model 
will  be  compared  to  results  of  experiments  on  rotational 
levels  of  two  typical  nuclei. 

The  predictions  of  this  model  for  the  nuclear 
radius,  the  moment  of  inertia  and  the  quadrupole  moment 
will  also  be  compared  to  experimental  values  of  these 
quantities . 

The  chapter  will  be  concluded  by  a  comparison  of 
our  model  to  some  others. 

The  Model 

The  model  we  want  to  quantize  will  consist  of  a 
spinless  mass-point  subject  to  the  fixed  potential 

V (r )  =  a (r4  -  2r2r2)  .  (3.1) 

o 

An  equivalent  picture  of  the  model  consists  of  two 
mass-points  moving  around  a  common  mass  center. 

A  quantization  axis  will  be  chosen  to  pass  through 
the  coordinate-origin  along  the  azimuthal  axis  of  the 
coordinate  system. 
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The  Energy  Eigenvalues 

The  energy  eigenvalues  E  of  this  model  are  given 

n  36 

by  the  one  particle  Schrodinger  equation  for  the  potential 


(3.1)  : 


tt"  4-  2y 

uo  +72 

*  ft 


(En^-a(r 


2r2r2)  -5-M1+AL)  u  =  0  (3.2) 

°  2yr  z 


l :  angular  quantum  number 
y :  mass  of  the  mass-point 
n:  total  quantum  number 

U^iR^r,  where  R  is  the  radial  wave  function. 


This  equation  (3.2)  cannot  be  solved  in  a  closed  form.  The 

eigenvalues  we  need  are  those  corresponding  to  angular 

excitation  but  minimal  radial  excitation,  which  means  n  =0 

r 

(the  radial  quantum  number)  and  n  =  £,  as  n—n  +SL,  The 
approximation  method  chosen  for  E  o  has  to  be  consistent 

JL 1  Xj 

with  these  conditions  on  the  quantum  numbers. 

Before  going  further  it  is  interesting  to  notice 

4  4/3 

that  in  the  fast-rotation  limit  when  V(r)  ar  ,  E  £-►...£ 
which  is  the  natural  limit  of  the  V.M.I.  model  (Moszkowski, 
79  p.  2394) .  Derivation  of  this  limit  for  E  is  in 
Appendix  B,  and  is  further  supported  by  derivations  of  E^ 
according  to  the  W.K.B.  method  (which  is  valid  for  n^  >>  £) 
in  Appendix  A. 

As  we  are  interested  in  finding  E  „  -  E  (where  E  is 

J  n£  o  o 

the  zero  point  energy  i.e.  n  =  £  =  0)  for  nr  =  0,  rather  than 

E  ,  we  can  use  the  following  approximation: 
n  36 


* 
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E  -  E  *  g(r  )  =  a(r4  -  2r2r2  )  +tL_JdA±ll 
n£  o  J  o£  oi  o  o V  ~  2 


(3.3) 


2yr 


Oj i 


where  r  is  defined  by  the  following  minimum  condition  on 


g  (r)  : 


=  0  .  (3.4) 

g(r)  is  here  the  effective  potential,  that  is  the  potential 

V(r)  (3.1)  with  the  centrifugal  term  added. 

The  centrifugal  term  is  not  modified  here  by 

2 

replacing  £(£+1)  by  (£+%)  ,  as  we  want  to  have: 


_d_ 

dr 


g  (r) 


r 

o 


r 

Oil 


£=0 


This  approximation  method  works  well  for  the  har- 

2 

monic  oscillator  when  (£+%)  '  dependence  is  assumed,  where 
it  gives: 


E 


£ 


{i+h)  instead  of  the  exact  result: 


E 


£ 


E  *-...£. 

o£ 


So  we  can  expect  the  method  to  give  close  to  correct 
eigenvalues  for  high  £. 

The  minimal  condition  (3.4)  is  equivalent  to  finding 
the  real  root  of: 


3 

x 


2  2 
r  x 
o 


fi2£(£+l) 

4ap 


0  . 


where : 


(3.5) 


x 
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The  Moment  of  Inertia 

The  moment  of  inertia  I  is  defined  as: 


I  = 


2  _ 

r  dy 


where : 


dy  =  dvp (r, ft) 


(3.6) 


(3.7) 


dv:  volume  element 

dy:  mass  element 

p(r#ft):  mass  density  or  probability  distribution 
of  the  nucleus 
Zl  :  angular  variables 


As  the  model  is  axially  symmetric  we  have: 

,  |2  ,U£,2 

P  £  ( r ,  }  =  li[Yj_0l  |-j!  (3.8) 

Y^  :  spherical  harmonics 
y  :  the  total  mass 


and  therefore: 


Xl  =  y 


U 


Z 


2  2 
r  dr 


Y*  sin^G  Y 0  d^ 
Zo  Zo 


The  radial  integral  that  yields  the  expectation  value  of 

2  2  2 

rQ:  <r0>  will  be  approximated  by  r  .  By  doing  the  angu- 

36  36  O  36 

lar  integration  for  Z  =  2  we  get: 


0.48 


(3.9) 


The  Quadrupole  Moment 


The  quadrupole  moment  of  the  model,  as  it  is  rotat- 
ing  around  a  symmetry  axis,  can  be  defined  as 


Q  = 


zp  (r,ft)  r2  (3  cos20  -  1)  dv 


(3.10) 


z:  number  of  unit  charges  e  assigned  to  the  mass- 
point 


where : 


z p ^  (r,  ft ) 


u 


=  z  Y 


% 


lo 


(3.11) 


is  the  probability  distribution  of  the  charge  ze  in  the 

2 

model.  By  making  the  same  approximation  for  <r  >  as 
before  we  get  for  1  =  2  the  following  quadrupole  moment 


~  2  4 

°2  Zro£  7  * 


(3.12) 


Comparison  with  Experiments 

Now  we  want  to  see  if  the  rotational  bands  of  our 

model  will  have  the  same  structure  as  the  rotational  bands 

170  238 

of  two  typical  nuclei  Hf  and  U  .  The  V.M.I.  model 
does  give  a  good  fit  to  experimental  results  for  these  two 
nuclei . 

In  our  model  we  will  assume  that  all  the  nuclear 
mass  is  located  in  the  mass-point,  though  experiments  do 
indicate  that  not  all  the  nuclear  mass  does  participate  in 
the  rotation. 

By  doing  this  we  will  get  the  minimum  values  for  1^, 

Q_  and  the  nuclear  radius  r  „ . 

2  oi 


19 


Table  1 


Computed  and  Experimental  Values  of  E0  , 

°2'  V  ro£  f°r  Hfl/° 


i 

E 

exp 

E 

comp 

ro£ 

keV 

keV 

fm 

„  ,  c  KeV 

a  =  11.6  — T 

0 

0 

0 

2.21 

fm4 

2 

101 

141 

3.2 

r  =  2 . 2 1  fm 

4 

321 

355 

3.3 

o 

6 

641 

665 

3.4 

2 

t?  “  122.4  keV 

8 

1042 

1051 

3.5 

2y 

10 

1504 

1500 

3.6 

12 

2015 

2000 

3.7 

14 

2565 

2546 

3.8 

16 

3150 

3130 

3.9 

18 

3764 

3751 

4.0 

20 

4417 

4403 

4.1 

±  i 

^ 2  2  comp 

2  j 

^2  2  exp 

^2  exp 

2  comp 

MeV""1 

Me  V" 1 

barn 

barn 

40 

59.8 

^5-7 

4.2 

Sources:  E^,  1^  :  (Sayer  et.al,,  75). 

Q 2  :  (Scharf f “Goldhaber  et.al., 


75  p.  279) 


Table  2 


Computed 

and  Experimental 

Values  of 

238 

V 

Q2'  ±2'  ro£  for 

•  u 

E 

E 

r  „ 

£ 

exp 

comp 

o£ 

keV 

keV 

fm 

ke  V 

0 

0 

0 

3.00 

a  =  2.45  ^jr 

fm4 

2 

45 

43 

4.4 

r  =  3.00  fm 

4 

148 

161 

4.4 

o 

~  2  0 

rV\  2 

6 

302 

334 

t:  -  87.4  keV  fm 

2y 

8 

518 

553 

— 

Zero  point  energy 

10 

775 

810 

of  42  keV  was 

12 

1076 

1098 

subtracted  from 

14 

1414 

1414 

5.0 

E 

comp 

2 

2 

I0 

2  comp 

-o  In  Q. 

^2  2  exp  2  exp 

^2  comp 

-1 

-1 

MeV 

MeV 

barn 

barn 

105 

133.6 

^  11 

10.3 

Sources:  E^, 

1^  :  (Sayer  et.al.,  75). 

°2 

:  CScharff 

-Goldhaber  et.al 

.  ,  75  p.  279) . 

75  p.  279) 
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The  hamiltonian  of  our  model  is  invariant  with 
respect  to  rotation  of  180°  around  the  symmetry  axis,  so 
the  angular  quantum  number  £  will  occur  in  sequences  of 
even  numbers  (Bohr  and  Mottelson,  75  p.  8) . 

In  the  preceding  data  fit  we  have  used  the  follow¬ 
ing  constants: 

2 

y  =938  MeV/c  for  one  nucleon 
■fi  =  6 .58  x  10~22  MeV 
C  =  3.00  x 1023  fm 


This  fit  to  experimental  results  is  by  no  means  a 
"best  fit”,  but  is  rather  intended  to  show  that  our  model 
can  reproduce  the  general  trends  of  rotational  bands  of 
some  typical  nuclei.  A  "best  fit"  method  would  be  compli¬ 
cated  and  slow,  due  to  the  fact  that  the  Schrodinger 
equation  (3.2)  (a  boundary  value  problem)  does  not  have 
a  solution  in  a  closed  form,  and  approximation  methods 
for  boundary  value  problems  become  time  consuming  when 
they  have  to  be  iterated  many  times. 

As  expected  the  values  for  I  and  Q  are  lower  than 
the  experimental  ones  due  to  the  assumption  that  all 
nuclear  mass  participates  in  the  rotation. 

If  we  assume  that  the  nuclear  radius  for  the  ground 
state  (£  =  0)  can  be  calculated  according  to  the  empirical 


formula : 


r  =  R  A 
o 


1/3 


(3 . 13 ) 
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R:  ^1.2  fm 

A:  number  of  nucleons 

then  we  get: 

r  (Hf 170 )  ^  6.7  fm 
o 

TOO 

r  (U  )  v  7.4  fm 
o 

which  is  nearly  twice  the  radii  our  model  gives  with  the 
conditions  we  did  assume. 


Comparison  with  Other  Models 

Prakash  et.al.  were  able  to  get  a  very  good  fit  to 
experimental  energy  values  by  a  model  that  has  the  follow¬ 
ing  Schrodinger  equation  (Prakash  et.al.,  73) 


,  fi2  d2  ,  h2£(£+l)  k  {  ,2  ...  .  A  .  .. 

(“2“”  - 2  - - 2"'"  +  2  (p  “  p0)  -E)^(p)=0  (3.14) 


dp 


2pp‘ 


where  y  =  M/4  is  the  reduced  mass  of  a  system  of  two  mass- 
points  . 

This  model  is  a  quantized  version  of  a  model  by 
Gupta  (Gupta,  69),  where  the  nucleus  is  represented  by  a 
prolate  ellipsoid,  rotating  around  an  axis  perpendicular 
to  the  symmetry  axis.  The  semi-major  axis  has  the 
length: 

Ro (1  -  ]w  <3-15) 

where  Rq  is  the  average  radius  and  3  is  a  parameter  measur¬ 
ing  the  deformation  from  a  spherical  form.  If  the  ellip¬ 
soid  is  divided  into  two  halves  along  a  plane  perpendicular 
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to  the  symmetry  axis,  then  it  is  easily  shown  that  the 
distance  p  between  the  two  mass  centers  will  be: 


p  =  7-  R  (1  - 
K  4  o 


4tt 


3) 


(3.16) 


By  a  change  of  a  variable  the  Schrodinger  equation 
(3.14)  can  be  written: 


2  2 

(_  jL  + 

2BT  d32 


v-  2  C  7 

~~  £(£+1)  +  (3  -  30)  ~  E ) ip 


=  0 


(3.17) 


with 


3MR 

R  =  15  _ o 

T  32  8  7T 


Cm  = 


45  p2 
6  4  7T  o 


I„  = 


45  2MR0 


T  128 


(1  + 


4  7T 


3) 


(3.18) 


— =  —-=&) 
bt  V 


and  finally  from  (3.14)  we  have: 


E  0 - y  £  (£  +  1  )  +  g-  (pT-p  )  +  (n  +  — Jll  ( - g-g: - +  77) 

r'  2y p  y  p 


£ 


(3.19) 


with 


p  —  p^  +  £  ( £+1 ) 


*-2 

TL 


yKp 


3  ' 


In  their  fit  to  experiment  n^_  =  0  is  required,  and  the  zero 
point  energy  is  subtracted  from  ^ . 

Due  to  the  complicated  angular  dependence  of  the 
nucleus  in  the  Prakash  model,  Iq  (£  -  0)  is  calculated 


2  ‘ 


instead  of  I 
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Table  3 


Our  Model 


Comparison  of  ours  and  the  Praka sh  Model 
170 

for  Hf _ with  Experimental  Values 

Z  x°  z  xj_  i°2 

MeV-1  MeV"1  fm 

40  3.2 


Model  of  Prakash  63.4 
Experimental 


2.84 

59.8  6.6 


r  ^  is  here  half  the  distance  between  the  mass- 
02 

points  in  the  Prakash  model,  which  is  equivalent  to  the 
radius  in  our  model. 

The  Prakash  model  gives  a  better  result  for  I,  as 

in  this  model  the  ground  state  is  assumed  to  be  deformed, 

which  is  experimentally  verified  to  be  true,  but  a  strange 

condition  in  a  central  potential. 

The  nuclear  radius  in  the  model  of  Prakash  R  is 

o 

given  by  (3.16)  and  is  larger  than  half  the  distance 
between  the  mass-points  due  to  the  spheroidal  ground  state. 

By  assuming  that  not  all  the  nuclear  mass  is  rotat¬ 
ing,  and  r  ^  in  our  model  will  increase  in  value. 

3  2  o2 

The  main  difference  between  our  model  and  the 


Prakash  one  is  that  Prakash  and  coworkers  take  a  successful 
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model  developed  by  Gupta  CGaupta,  69)  where  a  harmonic 
restoring  force  is  pre-assumed  and  quantize  it.  In  our 
case  we  begin  with  two  equivalent  models,  the  fourth 
order  cranking-  and  the  V.M.I.  model,  that  give  a  good 
description  of  the  rotational  energy  spectrum.  Then  we 
arrive  at  a  potential  that  classically  does  give  the  same 
results,  and  later  quantize  the  model.  Our  method  has 
therefore  strong  similarities  to  so-called  "inverse 
methods"  in  physics. 

In  this  work  the  main  emphasis  has  been  on  the 
model  itself  rather  than  using  its  results,  as  we  could 
not  find  similar  approaches  to  the  problem  in  the  physics 
literature . 


CHAPTER  4 


SIMPLE  FIELD  THEORETICAL  MODEL  OF  NUCLEAR  ROTATIONS 


In  this  chapter  we  take  a  simple  second  quantized 
model  describing  particles  rotating  on  a  circle  of  a  fixed 
radius  (Bassichis  and  Foldy,  64)  and  compare  it  with  the 
Harris  model  (Harris,  65)  and  an  ingenious  classical 
model  by  Thieberger  (Thieberger,  70)  that  has  the  same 
expressions  for  I  and  E  as  the  Harris  model. 

After  this  comparison  we  introduce  a  more  general 
interaction  into  the  model  of  Bassichis  and  Foldy  so  its 
energy  eigenvalues  will  have  the  same  properties  as  the 
rotational  bands  of  typical  even-even  nuclei. 

This  chapter  is  concluded  by  observation  of  several 
symmetries  of  this  new  model. 


Description  of  the  Models 

In  the  model  of  Bassichis  and  Foldy  (Bassichis  and 
Foldy,  64)  there  are  N  bosons  of  equal  mass  and  zero  spin 
with  a  two  body  interaction  present.  Its  hamiltonian  is: 


(4.1) 


The  radius  is  fixed  at  a  value  R.  If  the  potential 


is  assumed  to  depend  only  on  | 0^  -  0j|  then  H  will  have  the 
following  form  in  the  second  quantized  version: 


26 
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H  = 


^2^.2  •  ,  , 

2  akak  ^  V^kl~k4^k  +k  k  4-k  ak  ak  ak  ak 

k  2mIT  k  k  kxk  k  k  1  4  kl  k2'  3  4  kl  k2  k3  k4 


(4.2) 


f 


where  a^  and  are  creation  and  destruction  operators  of 
angular  momentum.  By  assuming  that  all  but  the  two  lowest 
single  particle  excitations  are  negligible  and  that  v(0)=0 
and  v(l)  =  g,  then  H  takes  the  form: 


H  = 


i  t  t  i  i  2ft  t  ^ 

h0)(ana  +aoa  )  +  g(a.a  (a,a.+a  a  )  +a  a.a  +a'  a,  a_  ) 
1133  22  1133  213  2  13 


(4.3) 


the  state  of  no  angular  momentum  is  labelled  by  2,  and  the 
states  of  angular  momentum  of  magnitude  ±fr  are  labelled  by 
1  and  3 . 

This  system  has  two  important  constants  of  motion 
(Bassichis  and  Foldy,  64  p.  A936),  the  total  number  of 


particles  N: 


f 


N  =  Y  a  .  a 

.  **,  i  i 
i=l 


(4.4) 


and  the  total  angular  momentum  L  of  the  system: 


T  f  f 

-  a]3].  -  a3a3 


(4.5) 


If  the  number  of  particles  in  level  3  is  given  by  n, 


then  the  state  function  is: 


%(N-£) 


+  ■*  l 


n=0 


CnjN,£,n> 


(4.6) 


where  the  C  1  s  must  satisfy  the  normalization  condition: 
n  J 
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%(N-£)  ? 

l  I  C  l  =  1  .  (4.7) 

n=0 

Before  going  any  further  let  us  now  look  at  another  model. 

Thieberger  (Thieberger,  70)  has  invented  a  classi¬ 
cal  mechanical  model  to  describe  the  successful  expressions 
of  the  Harris  model: 

E  =  ~  o2  (I  +  3Co)2 )  (2.1) 

2  o 

I  =  (I  +  2Co2)  .  (2.2) 

o  . 


Ro 


Figure  1:  The  Thieberger  Model 

The  model  consists  of  two  cylindrical  reservoirs, 
the  smaller  one  rotates  at  a  radius  with  the  angular 
velocity  o,  around  the  stationary  bigger  one.  The  liquid 
in  the  reservoirs  can  flow  freely  through  the  narrow  tube 
connecting  them.  The  model  is  stationed  in  a  homogeneous 
gravitational  field.  Due  to  this  field  and  the  rotation 
the  potential  energy  of  the  system  has  the  form: 
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V  =  (AM  +  BM2/2) 


where  M  is  the  mass  in  the  small  reservoir  when  oo 
o 


M 

M 


(4.8) 

=  0. 


Energy  Eigenvalues  of  the  Simple  Boson  Model 

It  is  now  possible  to  find  the  energy  eigenvalues 
of  the  hamiltonian  (4.3)  of  the  Bassichis  and  Foldy  model. 
The  method  for  doing  this  will  be  used  for  a  more  general 
hamiltonian  of  our  model,  and  to  estimate  the  order 
of  magnitude  of  the  parameters  involved. 

For  the  energy  eigenvalues  E  we  have  the  following 
expression : 


|  H  |  =  <xp  |  E  [  xp>  .  (4.9) 

By  using  the  state  function  (4.6)  expression  (4.9)  will  be 
equivalent  to: 

CE  =  C  [£  +  2n  +  g(N  -  £  -  2n)  (£  +  2n)  ] 
n  n  3 

+gCn-l  [(N~£“2n  +  2)(N“£“2n+l)  U  +  n)n]^ 

h 

+gC  , ,  [(N  -  £  -  2n-l)  (N  -  £  -  2n)  (£  +n+l)  (n  +  1)  ] 
n+l 

(4.10) 

which  is  an  exact  difference  equation  for  E,  that  has  to 
be  solved  numerically  in  most  cases. 

For  each  £  there  are  several  eigenvalues,  which 
correspond  to  vibrational  modes  in  the  nucleus.  So  for 
each  £  we  take  the  lowest  eigenvalue,  representing  a  pure 


. 
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rotational  excitation.  The  lowest  eigenvalue  can  be  found 
by  a  program  involving  the  Fortran  library  subroutine 
EQRT1S,  to  find  the  lowest  eigenvalue  of  a  tridiagonal 
matrix.  If  the  corresponding  eigenvector  is  needed  the 
subroutine  EQRT2S  can  be  used. 

Actually  we  need  the  eigenvalues  with  the  zero  point 
energy  E  (£=0)  subtracted  from,  to  compare  with  experimen¬ 
tal  results. 

Bassichis  and  Foldy  (Bassichis  and  Foldy,  64  p.A937) 
do  use  a  Bogoliubov  transformation  to  estimate  the  energy 
of  the  ground  state  of  the  hamiltonian  (4.3),  we  can  use 
this  transformation  to  check  the  numerical  solutions  for 
E  (4.10)  in  the  cases  of  low  i's. 

i  ■ 

In  the  approximation,  a^  and  a^  are  approximated 

P 

by  N2,  then  the  hamiltonian  (4.3)  will  have  the  form: 

ft  ft  ft 

H_  =  (a.a.  +a,aJfiu  +  gN(a1a  +a,a, +  a.a,  +  a  aj  .  (4.11) 

B  11  3  3  11  33  13  13 


This  hamiltonian  (4.11)  can  be  diagonalized  by  the  follow¬ 
ing  canonical  transformation: 


where 


t 

al  =  Cal ~  Sa3 


f 

a^  =  Ca^  -  Sa^ 


(4.12) 


t 

a-^  =  Ca^  -  Sa^ 


f 

a^  =  Ca^  ~  Sa^ 


C  =  [(.ho  +  gN  +  e)/2e] 

S  =  [ (ho  +  gN  -  e)/2e] 

e  =  [1  +  ho 

ho 


(4.13) 
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so : 


Then  the  hamiltonian  (4.11)  takes  the  form: 


ho  ( 1  +  -  fro  -  gN 


+  *“(1  +  Mi,  (atai  +  a+a3. 


(4.14) 


So  for  low  £  the  eigenvalues  are  approximated  by: 

P  P 

E_  U)  =  fro(l  +  -^)  2  -  ho  -  gN  tfrb  (1  +  ?^)  2£  (4.15) 

B  fro  ^  no 

and  therefore: 

P 

EBU)  -  EB(0)  =  A.fi<0(l  +  .  (4.16) 

This  expression  is  linear  in  £,  which  shows  the  need  for  a 

more  general  interaction,  if  the  model  is  going  to  be  used 

to  fit  values  from  experiments  for  E. 

The  computed  eigenvalues  of  the  hamiltonian  (4.3) 

are  shown  in  comparison  with  the  ones  from  the  Bogoliubov 
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approximation,  and  experimental  values  for  Hf  ,  m 
table  4,  just  like  the  Bogoliubov  approximation  predicted 
the  eigenvalues  are  almost  linear  in  £. 


Table  4 


Computed  and  Approximated  Eigenvalues 

of  the  Simple  Boson  Model  in  Comparison 

170 

with  Experimental  Results  for  Hf 


N  =  170  g  =  0.8  keV  =  10  keV 


E 

comp 

eb 

E 

exp 

keV 

keV 

keV 

106 

106 

101 

209 

210 

321 

309 

313 

642 

405 

415 

1041 

8 


• 
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Comments  of  Bohr  and  Mottelson  on  the  Harris  Model 

Bohr  and  Mottelson  (Bohr  and  Mottelson,  75  p.  25 
Vol.  II)  note  that  expansions  for  the  energy  E  and  the 
moment  of  inertia  I,  in  terms  of  the  angular  velocity  co, 
occur  often  naturally  in  the  analysis  of  rotational  pro¬ 
perties  in  terms  of  the  response  of  nucleonic  motion  to 
rotating  potentials,  and  give  as  a  reference  (Harris,  65) . 

For  a  system  rotating  around  its  symmetry  axis  the 
expansions  can  have  the  following  form: 


2  4  6 

I  =  a  +  $co  +  yco  +  6  oo  +  .  .  . 


(4.17) 


(4.18) 


The  expansion  for  E  (co)  has  been  used  successfully  to  fit 
results  from  experiments  with  high  degree  of  accuracy 
(Harris ,  65 ) . 


co  is  proportional  to  the  angular  momentum  Z,  so  the 

2 

hamiltonian  actually  can  be  written  in  terms  of  £  : 


H 


rot 


2 

V...&  +  .. 


2 

If  the  value  of  co  associated  with  the  quantum  state  £  is 


obtained  from: 


2  2  2 
fi  £(£+!)=  co  I 


(4.19) 


then : 


H  .  =  £  (£+1 ) 

rot  2a 


UU+l)  )2  +  ..  . . 


(4.20) 


(Bohr  and  Mottelson,  75  p.  25). 
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From  this  one-body  hamiltonian  (4.20)  and  the 
Thieberger  model  we  can  deduce  what  kind  of  interaction 
could  be  added  to  the  many-body  hamiltonian  (4.3)  for 
trial . 


Model  for  Nuclear  Rotations 

As  mentioned  before  Bassichis  and  Foldy  (Bassichis 
and  Foldy,  64)  use  the  hamiltonian: 


N 

H-  l 


h‘ 


i=l  2mR 


2 

ser 


+ 1  i 


v( 


e . 


e .  I) 


(4.1) 


which  describes  N  bosons,  on  a  circle  of  radius  R,  with  a 
two  body  interaction  V,  to  derive  the  second  quantized 
hamiltonian : 

ft  11*1*  ft  4"  2 

H  =fia)  (a^a^  +  a^a^)  +  g  (a2a2  (a^a^  +  a^a^ )  +  a^^a^  +  a2  a^a^)  • 

(4.3) 

The  classical  hamiltonian  for  particles  on  a  circle  inter¬ 
acting  through  the  two  body  interaction  V  is 

N  p0 .  , 

H  =  I  - i  +  j  l  v  (  1 0.  -  0  .  |  )  .  (4.21) 

i=l  2mR  i^j  3 

With  the  Harris  and  the  Thieberger  models  in  mind  we  want 
to  generalize  the  interaction  in  (4.21),  by  including  an 
interaction  that  does  depend  on  the  angular  momentum  of  the 
interacting  particles. 

If  this  interaction  is  an  even  function  of  the 
angular  momentum  then  the  resulting  hamiltonian  can  be 
expanded  in  terms  of  Pq . 


as  follows: 
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N 

H=  I 


P? 


- 2  +  I  J.  V'l0i  -  ei  I  >  +Y  2  l 

i=l  2mR  i^j  ij 


pe.pe . 


2mR' 


Y$  l 

i  jk& 


p  p  p  p 
0.0  .  0,  0O 
_ k  *■ 

„  2„4 
4  m  R 


+ 


(4.22) 


3 

In  the  first  quantization  PQ  is  replaced  bv  -  ih  ttz- /  this 
operator  is  hermitean  due  to  the  periodicity  of  the  wave 
functions.  A  product  of  two  hermitean  operators  is  also 
hermitean  if  they  commute.  We  have: 


3  01  '  Jtl  3  0  2  ^  0 

so  the  total  H  (4.22)  is  hermitean. 


(4.23) 

(See:  Messiah,  59, 


p.  256) . 

The  second  quantization  is  now  performed  by  intro- 

4. 

ducing  the  field  operators  ip  and  ^  1  ,  where: 


k 


1 

/2  7T  R 


-ik0 
e  a, 
k 


(4.24) 


The  first  two  terms  of  the  hamiltonian  are  derived  by 
Bassichis  and  Foldy  (Bassichis  and  Foldy,  64)  by  assuming 
that  v(k^  -  k^ )  =0  for  all  cases  except  when  k-^  ~  k^  =  ±p* 
For  the  two  body  interaction: 


Y 


I 

•  * 

II 


P  P 
0.0. 
1  J 

2mR^ 


the  matrix  element  is: 


1 


'fi2  h  ^k3\4  I  1  *k142> 


2m  E? 


1  fi2  ..  klk2 

2  2  Y  2 
2mR  (2rr) 


dQ^dO^  e 


i (k2  k3) e1  i  0^  k4) e2 


t  f 

3.1  9i  cl.,  3., 

kj  k2  k3  k4 


klk  2  2 

.  p2  ^  Y  °k~,k9  ^k,  ,k  a'  a}  a  a 
4mR  3  2  1  4  k]_  k2  k3  k4 


If  1/- 

12  *2  +  + 

- -j  *h  y  a,  a,  a.  a, 

4mR  K1  K2  k2  K1 


(4.25) 


for  the  three  level  system  this  will  be  (as  k.  =  ±k,0): 


k2h2  /  t  f  0  t  t  ,  t  t  s 

- 7y  y  (a  a  a  a  ~  za..a,s„a3  +  a-a»a„a.)  . 

A  TV  ^  JL  JL  ±  -L  JL  J-  d  s3  d  J  J  d 

2mR 


(4.26) 


The  matrix  element  for  the  pair  interaction  is: 

1  0  H  ,  ,  ,  ,  |  S  3  3  Si,,  ,  , 

"4  Y$  2  4  'K8K7k6k5  ^  807  77T  767  867  ^  k4k3k2k  1' 

4m  R  12  3  4 

(4.27) 

and  in  much  the  same  way  as  before  this  will  be  (as  we  get 

6  ,  6,  ,  <5V  ,  6,  ,  )  for  the  three  level  system 

k8'  4'  k6'  3'  6'K2'  5'  1 

.  4  4 

-  y3  ~~Y~2  (nSn43H-6N2N2  -4(.N2N  +NlN2)).  (4.28) 

4m  RJ  J 

If  we  set  now: 

,2,  2  _  -4,  4 

„  k  fi  ,  F  0nk  ,  .  . 

gF  =  - —  y  and  g-=  y0  — (4.29) 

2mR  4m  R 

then  the  added  interaction  can  be  written  more  simply  in 
terms  of  the  total  angular  momentum  i  as: 
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f(Jl)  =  gF(P  -  |)t4) 


(4.30) 


The  total  second  quantized  hamiltonian  of  our  model  will 
then  be: 


H  =  fto)CN1  +  N3)  +  g(N2(N1  +  N3) 


(4.31) 


where : 


f  NX  N3* 

This  interaction  has  to  be  an  even  function  of  £  in 

order  to  conserve  several  symmetry  properties  of  the 

hamiltonian,  as  will  be  shown  at  the  end  of  this  chapter. 

The  interaction  f(.£)  does  not  change  the  constants 

of  motion,  as  £  is  one  of  them  (see:  Dirac,  78  p.  92). 

In  table  5  the  energy  eigenvalues  of  H  (4.31)  are 
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compared  to  experimental  results  for  Hf  ,  and  m  table  6 
2  3  8 

for  U  ,  with  some  given  choice  of  the  free  parameters 
-hco,  g,  F,  a. 


. 


Table  5 
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Enerty  Fit  for  Hf 


N  =  170 

F  =  18.2 

in 

• 

i — 1 

II 

3 

keV 

a  =  1300 

g  =  0.8 

keV 

% 

E 

exp 

E 

comp 

keV 

keV 

0 

0 

0 

2 

101 

101 

4 

321 

314 

6 

642 

630 

8 

1042 

1040 

10 

1504 

1528 

12 

2015 

2076 

14 

2565 

2658 

16 

3150 

3249 

18 

3764 

3814 

20 

4417 

4318 

The  eigenvalues  E  of  this  model  are  within  4% 

3  comp 

of  the  experimental  values,  and  this  fit  is  by  no  means  a 
"best  fit"  of  any  kind. 


. 

Table  6 


Energy  Fit  for  U 


238 


N  =  238 

F  =  19.4 

fio)  =  0.1 

keV 

a  =  1550 

g  =0.4 

keV 

£ 

E 

E 

exp 

comp 

keV 

keV 

0 

0 

0 

2 

44.9 

44.9 

4 

148.4 

14  9 

6 

307.2 

310 

8 

517.7 

522 

10 

775.4 

779 

12 

1076 

1073 

14 

1414 

1391 

This  fit  is  within  1.6%  of 

the  experimental 

and  still  is  not  a  "best  fit". 

To  get  a  better 

fit  for  Hf 

a  higher  orde 

in  £  may  be  needed  for 

f  (£)  . 

I 


. 

< 

Moment  of  Inertia 


In  this  model  the  radius  R,  at  which  the  particles 
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are  rotating,  is  a  free  parameter.  It  can  be  fixed  for 
example  according  to  the  empirical  formula: 


R  =  R  A1</3 
o 

(4.32) 

R  -  1.2  fm 


The  state  function  of  the  model  (4.6)  is: 


Js(N-a) 

<P  =  l  c  jN,JL,n>  (4.6) 

n=0 


so  the  moment  of  inertia  I  will  be: 


I 


9  *s(N-A) 
mR2  l 
n=0 


n 


(4.33) 


where  m  is  the  mass  of  an  individual  nucleon. 

From  the  normalized  eigenvector  for  the  fit  of 
Hf  we  get: 


Mn-&) 

i 

n=0 


n  ^  39 


for  i  =  2 


(4.34) 


so  with  R  according  to  (4.32)  we  have: 

72  1l=2  r  83  MeV_1 

T1 

while  the  experimental  value  is  ^60  MeV  3 ,  which  we  can 
get  by  assuming  that  the  radius  is  5.7  fm  instead  of 


6.65  fm. 
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Interpretation  of  f(ft) 

It  might  be  interesting  to  see  what  kind  of  an 

✓ 

interaction  f (ft)  is  in  terms  of  the  original  operators 


Ni  =  a!ai 


N3  ^3^3 


We  know  that  ft  =  so: 


ft2  -  N2  -  2NxN3  +  N2 


d  =  N4  -  4N^N3  +  6N2N2  -  4NlN3  +  N4 


The  interaction  f ( ft )  can  then  be  written  as: 


f  U)  =  gFU2  -  ~  S.4)  =  gF(N2  +  n2  1 


n7  n. 


a  a 


2N1N3 


-  N2N2  +  -  (N2N,  +  NN2)) 
a  13  a  13  13 


(4.35) 


We  see  now  that  the  interaction  can  be  divided  into  a  two 
body-  and  a  pair  interaction  both  depending  on  angular 
momentum.  This  can  be  pictured  in  the  following  way,  where 
a  circle  signifies  a  particle  and  an  arrow  the  angular 
momentum  vectors. 


op  .  00 

repulsive 


9° 

attractive 


Figure  2: 


The  Two  Body  Interaction 


So  the  two  body  interaction  has  the  form  of  a  pairing 
force  to  decrease  the  total  angular  momentum. 


0000 


t 

oooo 


, ,  coco 

OOOO  '  -| 

*  *  oooo 


>  ^  * 


*4'  A/ 

attractive 


repulsive 


> ' 


Figure  3:  The  Pair  Interaction 


So  the  weak  pair  interaction  breaks  up  pairs  to  increase 
the  angular  momentum  of  the  model. 

On  the  level  of  the  original  operators  the  expan¬ 
sion  in  Z  has  a  reassemblance  to  a  multi-polar  expansion 
of  a  force  field. 

Symmetry  Properties  of  H 

Classically  time  reversal  is  defined  by: 

t  ->  -t  x  ->  x 

p  ->•  -p  L  **  ~L 

P  is  the  linear  momentum  and  L  is  the  angular  one 
(Frauenfelder  and  Henley,  75  p.  418) . 
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In  the  classical  hamiltonian  (4.22)  L  ,  which  is 
equivalent  to  P^,  occurs  only  in  even  powers,  so  these 
terms  are  not  affected  by  time  reversal.  We  can  therefore 

conclude  that  H  is  invariant  under  time  reversal . 


Parity 

The  conservation  of  parity  yields  a  quantum  number 
with  no  classical  analogue  (Frauenfelder  and  Henley,  75 
p.  359),  so  we  take  a  look  at  H  (4.28)  in  the  first  quan¬ 
tized  version. 

The  parity  operator  reverses  the  sign  of  a  polar 
vector,  but  leaves  an  axial  one  unchanged: 

x  ->  -x 

P  ->  -p 

L  +  L  so  P0  P0 

So  the  hamiltonian  H  will  preserve  the  parity  of  the  wave 
functions . 

Angular  Momentum 

We  know  that  the  angular  momentum  &  for  the  second 
quantized  hamiltonian  is  a  constant  of  the  motion,  and  we 
expect  it  is  also  for  the  more  general  classical  hamil¬ 
tonian  (4.22)  . 

If  L  is  conserved  for  H  .  then 

cl 


. 
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With: 


{Hcl'L} 


fHCl'L} 


=  0 


where  £  =  £ 

i  i 


I  ( 

k 


3H 


cl  3L 


3H  3L 


3  0, 


3P 


3Pq  30. 
6k  k 


) 


(4.36) 


we  have  that: 


3  £ 

3  0i. 


=  0 


3  £ 


3P, 


=  1 


(4.37) 


So  to  prove  that  the  Poisson  bracket  is  equal  to  zero,  we 
must  prove  that: 


3H 


cl 


3  0, 


=  0 


(4.38) 


which  is  equivalent  to: 


i  ik  v,ie‘"e)11 


=  0 


(4.39) 


As  the  potential  V  is  symmetric  in  y ®  j  aroun<3  the 
origin  we  can  expand  V(]yj)  in  a  cosine  Fourier  series: 


V 


(  y .  .  )  =  y  a  cos  ny .  .  . 
1  'll  1  L  n  'ii 

j  n  J 


(4.40) 


And: 


367  C°S  nYki 


367  cos  n^ki 


So : 


l  .1.  wr  v(l0i 

k  k 


0  .  |  ) 
J 


=  0 


and  therefore: 


(H  , ,L>  =  0  . 

cl 
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So  just  like  in  the  second  quantized  version  the  symmetry 
of  the  potential  ensures  the  conservation  of  the  angular 
momentum. 

The  potential  V  does  not  depend  on  velocity,  and 
Hcl  <^oes  not  depend  explicitly  on  t,  so  the  hamiltonian 
is  a  constant  of  the  motion  and  equals  the  total  energy 
E  (Goldstein,  77  p.  220). 
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Appendix  A 


E  ^  according  to  the  W.K.B.  Method  (nr>>£) 

In  this  appendix  we  want  to  derive  the  energy  eigen¬ 
values  of  the  Schrodinger  equation  for  the  potential  V(r)  = 
4 

ar  ,  by  the  W.K.B.  method  (Morse  and  Feshbach,  53  p.  1101). 
The  Schrodinger  equation  is: 


+  4  ut  =  ° 


with : 


qr  =  (e  -  v  (r) 


£  (£+1) 


) 


U£ .  R£  r 


where : 


Rn  :  the  radial  wavefunction 

£ 

£  =  2yE/h2 

v  ( r )  =  V  ( r ) 
fi 

]i  :  the  reduced  mass  of  the  system  . 


( A.  1 ) 


(A. 2) 


By  making  the  transformation  r  =  e 


x 


equation  (A.l)  takes  the  form: 


x"  +  e2x  (e  -  v  (eX)  -(£  +  -)  e  2x)  x  = 


and  U  =  e  2X  x  the 


0  .  (A.  3) 


This  equation  (A. 3)  can  be  treated  in  the  same  way, 
by  the  W.K.B.  method,  as  the  one  dimensional  Schrodinger 
equation,  because  the  singularities  are  now  at  x  =  +  00. 
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The  condition  for  a  bound  state  according  to  the 


W.K.B.  approximation  is: 


x. 


x. 


qx  dx  =  <n  +  2>  *  . 


(A. 4) 


If  n^  >>  i,  x^,x^  are  positive  roots  of  q  x  which  is  the 
coefficient  to  y  eq.  (A.  3),  and  n^  is  the  radial  quan¬ 
tum  number. 


By  going  back  to  the  variable  x  we  have  more 


simply: 


( o  + \  ^  h 

(s  -  v  (r)  -  )  dr 


(nr  +  2 )  tt 


(A. 5) 


4  2 

For  v(r)  =  ar  ,  q  will  be: 


2  4  (l+h) 

q  =  e  -  ar  - - 4~ 

^r  2 


(A. 6) 


2 

r^,r^  will  be  two  positive  roots  of  q^  corresponding  to 


1'  3 
the  c 


lassical  turning  points  of  the  effective  potential. 


22  6  2  .  -  N  2  n 
r  q^  =  ar  -  sr  +  (Z+h)  -  0 


(A. 7) 


By  changing  the  variable  r,  in  (A. 7),  to  x2  we  get: 

2 


X3  -  £.  X  +  AIM 

a  a 


0  . 


(A. 8) 


Now  there  are  three  possibilities: 


1. 


If  A  = 


tt+k)  4 

4  a2 


27a' 


>  0 


then  there  is  just  one  real  root. 
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2.  If  the  same  quantity  equals  zero,  then  there  are 
three  real  roots,  two  of  whom  are  equal. 

3.  If  the  same  quantity  is'  always  larger  than  zero, 
then  there  are  three  real  unequal  roots  (Selby,  75  p.104). 
The  3.  possibility  can  be  rewritten  as: 


£  > 


(l .  61 ( i+h) ) 


4/3 


a 


1/3 


which  we  will  see  later  is  always  satisfied. 

To  see  what  we  can  expect  from  the  W.K.B.  method, 
we  can  easily  do  the  integration  for  Z  =  0 : 

4/e/a 


0 


A  ' 2 

(e  -  ar  )  dr  =  (nr  +  tt 


(A. 9) 


3/4  gives  a  better  approximation  here  than  1/2. 


Change  the  variable:  r  =  ax,  x  =  —  r,  a  =  4/e/a  , 

Uv 


then  we  have 


(n  +  7-)  tt  = 
r  4 


which  comes  to: 


(e  -  ex^)  2  dx^  /e/a 


0 


(A. 10) 


e  = 


,  .  3  x  t  ^  4/3 

{nr  4}  71  a4 


1  a  % 

J  (1  -  x  )  dx 
0 


(A.  11) 


I  = 


1  j  1 

(1  -  x4)  2  dx  = 


1  -  x 


dx 


0 


0  A  - 


From  the  integral  tables  in  (Selby,  75)  we  have: 


I 


A—  k  (i)  -  0.875 

3/2 


where  K(k  )  is  the  complete  elliptic  integral  of  the  first 
kind . 

Expression  (A. 11)  can  now  be  written  in  the  final 

form: 


0  - 


/  ^3.  k]4/3 

(nr  +  4)  tt  a4 


0.875 


(A. lib) 


According  to  (Selby,  75  p.  104)  the  roots  of  q  can  be 
expressed  as 


If: 


x^  =  A  +  B 

X3  =  -  ~(A  +  B)  -  (A  -  B) 

X 2  =~i(A  +  B)  +  (A  -  B) 


f  \ 

A 

v _ 

\ 

a+h) 

2a 


±  ( 


U+h) 

4a^ 


3  t 


27a 


1/3 


(A. 12) 


(A. 13) 


Im  x 


Figure  4 :  Location  of  the  Roots 

of  <2, 


The  relative  lo¬ 
cation  of  the  roots  is 
shown  in  figure  4 . 

The  arrows  show 
the  trends  of  the  roots 

as  e  increases  com- 

2 

pared  with  (Jl+%)  . 


lx 


So  in  order  to  obtain  the  energy  eigenvalues,  we 


have  to  evaluate  the  integral: 
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I  = 


,  4  a+h)2,"2  , 

(e  -  ar  -  — — )  dr 


(A. 14) 


This  can  be  done  in  terms  of  elliptic  functions,  if  we 
make  the  following  changes:  first  the  integral  can  be 
written  as: 


I  = 


/  2  6  /rt  ,  i  v  2,  2  dr 

(er  -  ar  -  a+h)  )  — 


(A. 15) 


2  h  1  -h 

then  the  substitution:  r  =  x,  or  r  =  x  ,  with  dr  =  2-x  dx 


will  give: 


I  = 


1 


(ex  -  ax3  -  U+!s)V 


(A. 16) 


By  multiplying  and  dividing  the  integrand  by  the  square 
root  in  it,  and  taking  out  a  factor  of  l//a  we  get: 


I  = 


ex  -  ax 


a+h)2 


2  /a  J  /  e  3 

X  /—  x  -  x  - 
J  a 


dx 


(A. 17) 


a+h) 


The  polynomial  inside  the  square  root  has  the  root's 
X1 '  x2  and  x  .  In  order  to  write  I  in  the  standard  form 
for  elliptic  integrals  we  have  to  write  the  integrand  in 
terms  of  the  roots  x^ ,  x2  and  x^ . 

The  limit  of  integration  is  between  the  positive 
roots  x^  and  x^  so  I  will  have  the  form: 


I  = 


2/a 


1 


x. 


ax 


dx  - 


x  ^  (xi~x)  (x“x2  ^x~x3^  2/a 


2  1 
a+h)  ~ 

-X. 


dx 


x  /  (x1~x)  (x-x3)(x-x3) 


= J1  +  I3 


(A. 18) 


. 


/  '  ■  ; 
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To  do  the  integration  in  I,  we  use  the  general 
formula  (Byrd  and  Friedmann,  71  p.  81). 


U. 


R(t)  dt 


^Ta-t)  (t-b)  (t-c) 


=  9 


R (a  -  (a-b)  Sn  U)  dU 


(A. 19) 


0 


which  is  valid  for  y  in  the  range  a  >  y  >  b  >  c,  R(t)  is  any 
polynom  in  t.  Other  symbols  are  defined  as  follows: 


g  = 


SnU. 


/a  - 


=  /  a  -  y 
J  a-b 


k  '  the  parameter  of  the  elliptic  functions  is  defined  as 


a-b 
a  -  c 


In  our  case  we  have 


=  a 


then  SnU^  =  1 


*3  =  b  =  y  z=>  U1  =  K (k  )  =  K 


x2  =  c 


where  K  is  the  complete  elliptic  integral  of  the  first 
kind. 


can  now  be  written  as: 


K 


I = 


^  ^1  X2  0 


(e  -  a  (x^  -  (x^-x^  )  Sn2U)  2  )  dU 


/a/Xf  -  x2 


( s  K  —  a 


K 


0 


(x2  -  2x1  (x1~x3  )  Sn2u  +  ( x ^-x  ) 2  Sn^U )dU 


/a/x^  -  x2 


(£K  -  I2)  . 


(A. 20) 
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From  (Byrd  and  Friedmann,  71  p.  191-192)  we  have: 


I2  =  a(Kx2  -  2x1(x1~x3)  ~  (K-E(K)) 

+  (x  -x  )2  -~  (  (2  +  k2)K  -  2  (1  +  k2)E(K)  )  (A. 21) 

1  J  3k 

as  E(0,k)  =  0,  SnO  =  0,  CnK  =  0.  E  (K)  is  the  complete 
elliptic  integral  of  the  second  kind.  Later  it  will  be 
shown  that  the  coefficient  to  E(K)  is  equal  to  zero. 

Now  we  have  to  solve  the  integral: 


I3=- 


2  A 


?1 


x. 


(£+%)  dx 


x 


/  (x,  -  x)  Zx  -  x0 )  (x  -  x_) 


In  our  method  to  solve  this  integral  x^  must  be  larger  than 
zero  by  some  small  amount  at  least. 

To  evaluate  I  we  integrate  along  the  complex 
contour  C  shown  in  figure  5. 


Im  x 


Figure  5:  The  Contour  C 
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The  integration  along  C,  I  will  yield  zero  as  no 

pole  is  inside  C.  The  integrand  vanishes  when  the  inte- 

* 

gration  in  the  upper  part  of  C  is  taken  to  infinity.  So 
we  have : 


x 


+  J' 


'D 


r  3 


+  I. 


ri  Res  f 


x=0 


(A. 22) 


ft  means  here  principal  value  and  Res  f|  indicates  that 

I  X“ u 

the  residue  value  will  be  taken  of  the  integrand  at  x  =  0. 


We  know  that  I~  will  be  real  valued  and  1=0,  so 
v  3  00  C 

r3 

as  the  integrals  J  and  are  purely  imaginary,  they 

*2  x 

must  cancel  each  other.  Then  we  have 


0  = 


+  I  -  ir  Res  f 


—  oo 


x=0 


x, 


—  00 


+  J3  +  f  U  + 


So  now  we  can  express  as: 

x^ 


I3- 


2/a 


U+h)  dx 


7T 


^  x/ (x^-x)  (x-x2)  (x-x^) 


U+h)-  (A. 23) 


To  evaluate  this  integral  we  have  from  (Byrd  and  Friedmann, 
71  p.  70)  the  general  expression 


U. 


I'  = 


R(t)  dt 


—  00 


/  (a-t)  (b-t)  (c-t) 


=  g 


R  ( a  ~  (a-c)  nS  U)  dU  .  (A.  24) 


0 


I 


R(t)  is  any  rational  function  of  t.  The  condition  for  y 
is: 

* 

a  >  b  >  c  >  y  . 

And : 

g  -  - - -  U  =  K 

/a  -  c 

in  our  case  we  identify  the  roots  as  follows: 


X1  =  a 


x2  =  b 


x3  =  c 


then  can  be  expressed  as: 


K 


t  -  ( &+h) 

3 


2/a  /x,- 


x 


dU 


1  3  0 


(xn  -  (x  -x_)  nS  U) 


|  U+h) 


(A. 25) 


1  2 


If  we  now  define  the  parameter  y  =  x^/(x^~x^)  which  is 
always  less  by  value  than  one,  then  the  integral  itself 


can  be  written: 


-y 


K 


0 


Sn  U  dU 
2  2 

1  -  y  Sn  U 


with 


nS  = 


Sn 


(A. 26) 


on  the  condition  that 


2  2 

k  <  y  <1  where 

This  integral  has  the  solution 


2  X1  "  X3 


Xl"X2 


K 


-  y 


Sn2UdU 

2  2 

1-y  Sn  U 


-  y 


o  irCl-A  (0,k)) 

z _ o _ 

2/y2 (1  -  y2 )  (y2  -  k^7 


(A. 27) 


0 
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with: 


A  .  -1  /  1  -  y2  .  -1  /  x2 

0  =  sm  / - * —  =  sm 


,2 


X2  “  X3 


AoC0,k)  is  an  elliptic  integral  of  the  third  kind  (Byrd  and 
Friedmann,  71) . 

After  a  simple  algebraic  manipulation  can  be 
written  as: 


TT 


I3  2  /- 


U+h)2  (1-V0'k)) 


TT 


/a/Xj-x^ 


/  X1X2X3 
Xl“X2 


-  ~  a+h) 


(A.  28) 


2  .  .  2 

In  the  polynom  q  the  coefficient  to  x  is  zero, 

X 

which  imposes  the  following  condition  on  the  roots: 


X1  +  X2  +  X3  =  0  • 


In  the  integral  the  coefficient  to  E(K)  is: 

2 

a2x  a  (x-,  -  x^ )  ? 

- (v  ~  x _  )  +  - - - ^ -  2  ( 1  +  k  ) 

k  1  J  3k 


which  can  be  proven  to  equal  zero,  due  to  the  sum  condi- 

2 

tion  on  the  roots,  and  the  fact  that  k  =  x^-x^/x-^-^ . 

The  sum  condition  enables  us  as  well  to  eliminate 


:3  from  all  our  expressions.  If  now  the  factor  to  K(kz) 


in  integral  I  is  simplified,  I  will  be: 


I  =  I1  +  I3 


/a  (x1-x2Y 


k(e+j  (x2  +  x^x2  -  x2)  } 


2  ( 1  -  A  ( 0 ,  k )  ) 


/ 


/Xlx2(x  1+X2) 


X1  “x2 


(A. 29) 
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And  the  condition  for  a  bound  state  in  the  potential  is, 
as  stated  before: 

* 

I  -  (nr  +  i-)  it  =  0 

Now  we  have  to  check  if  this  expression  for  bound 
state  energy  values,  when  we  require  that  £  -  0  will  yield 
the  same  eigenvalues  as  when  £  =  0  was  initially  assumed 
for  V 

ef  f 

If  £  (£+1)  dependence  had  been  assumed  instead  of 
(£+%) ,  then  we  would  have: 

x1  +  /e/a  x2  -*  ~/e7a 

x^  ->  0 +  as  £(£+1)  ->  0 

therefore:  ->  1/2 

and  the  bound  state  condition  becomes: 


/  i 


1 


A 


K  {e  +  |  (-  f) 
3  a 


-  (nr  +  J)w  =  0 


2/e/a 

which  when  rearranged  will  give: 


,  ,  3»  1/4 

nr  +  4"  3 


3  7T 


/ 2  K(h) 


4/3 


,  ,  3,  p  4/3 

r  (n  +  -t-)tt  a  4  ^ 

— £ — ± - 1 

0.875  j 


(A. 30) 


(A. 31) 


3  1 

if  (n  +  7-)  is  assumed  instead  of  (n  +  7-)  in  bound  state 
r  4  r  2 


condition . 
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For  >>  i  equation  (A. 29)  can  be  written 
explicitly  for  E ,  as  fE/u  and  the  term  including 

can  be  ignored. 


E 


3 

(2n^  +  £  +  — ) 


ir/2 . 3 
4K(%)> 


4/3 


*2  2  > 
/n  x 

(2ir>  “ 


1/3 


E 


(A 


A 

o 


.32) 


Appendix  B 


E  in  the  Limit  when  £  >  1  and  n  =  0 
n£  ,  r 


for  V(r)  =  ar 


The  energy  eigenvalues  E  ^  -  Eq  (the  zero  point 

energy  Eq  has  been  subtracted  here)  can  be  approximated 

4 

for  the  potential  V(r)  =  ar  by: 


/  \  4  ,  ft2  (£  +  h)  2 

g(r  )  =  ar  +  -1 — 

o  o  „  2 


(B .  1 ) 


2'iTr‘ 


o 


g:  is  the  effective  potential:  that  .is,  the 

modified  centrifugal  term  has  been  added 

to  the  potential  V(r). 

r  : is  the  value  of  the  radius  at  which  g(r) 
o  ^ 


takes  a  minimum  value. 


if  £  >  1  and  n  =0.  The  minimum  value  r  can  be  shown 

r  o 


to  be: 


r  = 
o 


2  2.  !/6 
ft  u  +  %)2 


4ay 


(B .  2 ) 


If  this  expression  for  rQ  is  now  used  in  (B.l)  we  get: 


E 


n  =0 
r 


4/3  f,2  2/3  1/3  1/3  1/2 

-  g(r  )  =  U+h)  a  (f-  +2  ) 

O  2]i  4 


r(4+>s)4'/3  1.890 


*2  2 

“J 


1/3 


(B .  3 ) 
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